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We investigate theoretically the phase diagram of an insoluble charged surfactant monolayer in 
contact with a semi-dilute polyelectrolyte solution (of opposite charge). The polyelectrolytes are 
assumed to have long-range and attractive (electrostatic) interaction with the surfactant molecules. 
In addition, we introduce a short-range (chemical) interaction which is either attractive or repulsive. 
The surfactant monolayer can have a lateral phase separation between dilute and condensed phases. 
Three different regimes of the coupled system are investigated depending on system parameters. 
A regime where the polyelectrolyte is depleted due to short range repulsion from the surface, and 
two adsorption regimes, one being dominated by electrostatics, whereas the other by short range 
chemical attraction (similar to neutral polymers) . When the polyelectrolyte is more attracted (or at 
least less repelled) by the surfactant molecules as compared with the bare water/air interface, it will 
shift upwards the surfactant critical temperature. For repulsive short-range interactions the effect 
is opposite. Finally, the addition of salt to the solution is found to increase the critical temperature 
for attractive surfaces, but does not show any significant effect for repulsive surfaces. 

PACS numbers: 82.35.Gh, 82.35.Rs, 61.41. +e 



I. INTRODUCTION 

Complexation between surfactants and macro- 
molecules has^enerated ^^r^t d eal of research, J 30 ^ 

and experimentally^!!. f^Ull l^j^ffllinla l^f. 

This research is motivated in part by the role surfactant- 
macromolecule complexes plays in biological sys- 
tems |g, Il0| . where the adsorption of a variety of 
macromolecules onto charged lipid cell membranes is a 
precursor to many important biological processes. 

The complexation of charged surfactants and macro- 
molecules, besides its relevance to biological physics, 
carries interesting questions about the interplay of 
short-range interactions between surfactants and macro- 
molecules, as compared with long-range electrostatic in- 
teractions. We consider a multicomponent monolayer 
where the possibility of a lateral phase demixing be- 
tween the various constitutes exists and is characterized 
by a demixing (miscibility) curve in the temperature- 
concentration plane. For very long polyelectrolyte (PE) 
chains, the electrostatic interactions have two opposite 
effects on the adsorption. Electrostatic attraction be- 
tween surfactants and PE chains increases the complexa- 
tion, while the electrostatic repulsion between the PE 
monomers tends to decrease the complexation. This 
interplay also stimulated the interest in the adsorp- 
tion of polyelectrolytes on oppositely charged surfaces 

m m m es m El IM Hi in in mis in . 

Several authors have addressed a closely related prob- 
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lem of neutral polymers adsorbing on a surfactant mono- 
layer [UlEl- It was shown that for polymer-surfactant 
systems the surfactant demixing occurs for higher tem- 
peratures as compared to a pure surfactant system. 
Other works dealt with another similar problem of glob- 
ular proteins (rigid macromolecules with no flexibil- 
ity) adsorbed on charged membranes 0, 0- A recent 
paper also deals with the adsorption of one flexi- 
ble polyelectrolyte chain on a surfactant monolayer us- 
ing Monte Carlo simulations. The possibility that the 
polyelectrolytes can desorb entirely from the charged 
surface when the solution has sufficiently strong ionic 
strength Ej| was not previously discussed. 

In this paper, we analyze the demixing of surfactants 
on the monolayer plane using numerical solutions of the 
mean- field equations. We find that the short-range in- 
teractions between the PE and the surfactants have a 
greater effect on the shape of the phase separation curve 
than the long-range electrostatics. Three regimes for the 
surfactant-polyelectrolyte complexes are found: (i) a de- 
pletion regime, where the polyelectrolytes do not adsorb 
at all onto the surfactant layer; (ii) an electrostatic ad- 
sorption regime, where the main attractive interaction 
between surfactants and polyelectrolytes is the electro- 
static attraction; and (iii) a chemical adsorption regime, 
where the main interaction between the polyelectrolytes 
and the surface is short-ranged (chemical) and attrac- 
tive. In the first two regimes, the surfactant demixing 
curve (binodal line) decreases in comparison to that of 
the pure surfactant monolayer, while in the third regime 
the demixing curve increases. We also find that for the 
electrostatic adsorption regime, the adsorbed amount for 
a system in the coexistence region is almost the same as a 
monophasic system with the same average concentration. 
In the chemical adsorption case the adsorbed amount in- 
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creases more substantially due to the phase demixing. 

The outline of this paper is as follows: we start by 
introducing the free-energy functional for a PE solution 
with salt in contact with a charged surfactant monolayer. 
In Sec. IIIII we present our numerical results and discus- 
sion. In Sec. II VI we present a simple scaling analysis for 
three limits: (i) PE depletion from a short-range repul- 
sive surface; (ii) PE electrostatic adsorption to a short- 
range repulsive surface; and (iii) the short-range attrac- 
tive surface case. Finally, we present some concluding 
remarks in Sec. [3 



II. FREE ENERGY FORMULATION 

Consider an aqueous solution of polyelectrolyte chains 
and salt ions. The bulk PE solution is in contact with a 
planar and rigid interface, on which resides an insoluble 
monolayer of charged surfactants. This monolayer can 
be thought of as a Langmuir monolayer at the air-water 
interface, or as the outer leaflet of a lipid bilayer mem- 
brane. The charged surfactants are assumed to be able 
to move freely on the interface, but can not dissolve into 
the solution (see Figure 1). 

The free energy of such a system is composed of three 
contributions. The monolayer free energy F sur f , the solu- 
tion free energy F so \ and the interaction free energy Fi nt . 
We use a standard lattice-gas formalism Q to derive the 
demixing free energy of the surfactant monolayer: 



Rnd = b- 2 [cine + (1 - c) In (1 - c) (1) 
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where b 2 is the close packing area of the surfactant head 
groups. The first two terms account for the surfactant 
entropy of mixing, where c is the surfactant area frac- 
tion. The third term represents a short-range interaction 
between the surfactants with v~ x as the interaction pa- 
rameter. The fourth term is electrostatic contribution 
of the charged surfactant layer, where Z is the valency 
of the surfactant head group (assumed henceforth to be 
negative) and £ s = e^(0) / 'ksT is the surface potential on 
the monolayer plane in units of fc^T. 

Without the fourth term, F sur { is the free energy com- 
monly used for neutral Langmuir monolayers. For sim- 
plicity, we take it to be symmetric under the exchange 
of c and 1 — c. The non-convexity of -F sur f implies that 
there is a demixing region with two phases: (i) a dense 
phase with high surfactant concentration, and (ii) a di- 
lute phase with low surfactant concentration. The phase 
demixing curve is determined by: 
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where v = Vb{c) is the binodal line. Below a critical value 
v c = 1/2 and for v{c) < Vb{c) the monolayer phase sep- 
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FIG. 1: A schematic view of the polyelectrolyte-surfactant system. 
The polyelectrolyte solution contains polyelectrolyte chains and mono- 
valent salt (co- and counterions). The negatively charged surfactants 
are present only on the air-water interface, and phase separate into 
dense and dilute lateral phases. 



arates into two coexisting phases. Inside the coexistence 
region, another line of interest is the spinodal line, sat- 
isfying F s " rf (c) = 0. The spinodal line, separating the 
unstable and metastable regions with respect to surfac- 
tant density fluctuations, is given by: 



2c(l-c). 



(3) 



Its maximum coincides with that of the binodal line at 
the critical point, c c = 0.5 and v c = 0.5. Therefore, 
examination of the spinodal line gives a good initial ap- 
proximation to the phase behavior. 

The free energy of the bulk PE and salt solution, F so i = 



F V o\ + F[ c 



F c i, is more complex, and is composed of 



three parts: (i) the polymer free energy F po i; (ii) the 
small ion entropy -Fiona! an d (iii) the electrostatic free 
energy F e i. The polymer free energy is |30| : 



Fpoi 


poo 

/ dx 
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where the first term accounts for the polymer chain flexi- 
bility with a the monomer size, 4> 2 the bulk concentration 
of monomers and r) 2 (x) = 4> 2 /<j> 2 the local monomer con- 
centration renormalizcd by its bulk value. The param- 
eter v cx is the polymer excluded volume (second virial) 
coefficient. The second term originates from the excluded 
volume interactions between the monomers as well as the 
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equilibrium condition of the PE solution in contact with 
a bulk reservoir. 

The change in small ion entropy with respect to the 
reservoir is: 



k B T 



E 



d.r 
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where c± (x) is the local concentration of the positive and 
negative ions and c b is their bulk concentration. 
The last term of F so i, the electrostatic energy, is: 



F e l 



(6) 



d.r 



where the dimcnsionless electrostatic potential (rescaled 
by the temperature) is £ = eip/ksT. The first three 
terms are the interaction of the electrostatic potential 
with the positive ion, negative ion and polyelectrolyte 
concentrations, respectively, where the parameter / is 
the fraction of charged monomers on the PE chain. The 
last term is the self-energy of the electric field, where Ib = 
e 2 /eksT is the Bjerrum length. For water with dielectric 
constant e = 80 at room temperature T — 300K, Ib is 
equal to about 7 A. 

The last contribution to the free energy, the non- 
electrostatic interaction free energy -Pint, can be written 
as 0: 



1, 
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In the above rj s = ?/(0) is the renormalized monomer con- 
centration at the surface, x — 0. The strength of the PE- 
surfactant interaction is defined by a phenomenological 
parameter e ps , which has usints of volume. For positive 
values of e ps , the short-range interaction of the PEs with 
the surfactant is more favorable than those between the 
PE and the bare interface, while for e ps < the opposite 
occurs. At concentration c — c* the surface is indiffer- 
ent ( "special transition" line) Q, Q ■ For e ps > and 
low surfactant concentrations c < c* , the short-range in- 
teractions between the surface and the PE are repulsive, 
while for higher surfactant concentrations c > c* these 
interactions are attractive. 

In order to find the equilibrium state of the PE- 
surfactant system we minimize the total free energy, 
F = F BUr f + F pol + F ions + F e i + Fint with respect to the 
local small ion densities c± , electrostatic potential £ and 
monomer order parameter tl This minimization yields 
the following equations |3(J, LiH |33] : 



da; 2 



= k 2 sinhC + k 2 m (e c - rf) + 
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where k 1 = (87r/sc sa it) L/ * is the Debye-Hiickel length 
scale for the screening of the electrostatic potential in 

the presence of added salt and fc^ 1 = (47t/_b0 2 /) is 
the corresponding length for the potential decay due to 
counterions. Note that the actual decay of the electro- 
static potential is determined by a combination of salt, 
counterions, and polymer screening effects. 

The solution of Eqs. JSJl and © requires four bound- 
ary conditions. Two of them arc chosen at the bulk, far 
from the surfactant monolayer, and two at the surfactant 
layer. At the bulk we set the electrostatic potential as 
C(x — > oo) = 0, while the monomer concentration is set 
equal to the bulk monomer concentration r)(x — » oo) = 1. 
The other two boundary conditions are obtained by in- 
tegrating Eqs. © and @ from to a small distance 
near the surface incorporating the surface interactions of 
Eqs. ©, ©: 



dC 

d.r 

drj 
dx 



AttIbcZ 



b 2 

3£ De 



x=0 



b 2 a? 



'-t](0) 



(10) 

(11) 



Equation l|10|) is the usual electrostatic boundary condi- 
tion for a given surface charge density |oj = cZ/b 2 , while 
Eq. is the Cahn - de Gennes boundary condition [3^ , 
which is often used for neutral polymers 29] . 

Equations l|5)l-i|llD are solved numerically using the re- 
laxation method |33, as was described previously [33| . 
For each surfactant concentration, we calculate the value 
of v for which F"(c) = 0, corresponding to the spinodal 
curve. Then the full binodal line is calculated numeri- 
cally. We find that despite the electrostatic interactions 
between the surfactants and the PE chains, the main 
characteristics of the spinodal and binodal lines depend 
on the short-range (non-electrostatic) PE-surfactant in- 
teractions. There are two main regimes for the demixing 
curve. For large and positive e ps (c — c*) values (short- 
range attractive surface) , the demixing curve is shown to 
be always higher than the pure monolayer one, Eq. 
For large negative e ps (c — c*) values (short-range repul- 
sive surface) , the curve is always lower than the pure one, 
Eq. ©. 

In order to incorporate the temperature dependence 
of the various parts, we note that v, the inverse Flory- 
Huggins parameter, is proportional to the temperature T. 
The chemical interaction strength parameter e ps is pro- 
portional to T _1 . The other parts of the free energy have 
a weaker dependence on T, which can be ignored in this 
simple model . Their contribution is neither purely en- 
thalpic nor purely entropic, and is less significant. There- 
fore, in order to simulate the change in temperature in 
our simple model, we follow the method shown previously 
in Ref. 3] and choose v ■ e DS = to be constant. 
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III. NUMERICAL RESULTS 
A. The Three PE-Surfactant Regimes 

We find three distinct regimes for the PE-surfactant 
system. These regimes are presented in Fig 2, and dis- 
cussed more thoroughly in Sec. IIVI For repulsive sur- 
faces, very low c and high c sa it, the PEs deplete from the 
interface. The high amount of salt screens the electro- 
static interactions between the constituents of the solu- 
tion, and the free energy can be written as a sum of three 
terms: an interface term, a neutral polymer term and a 
salt solution term. For higher c, short-range repulsive 
surfaces and low salt, the PEs adsorb electrostatically on 
the surface. In this case, the free energy is dominated by 
the electrostatic attraction between the PE chains and 
the surfactant layer. For the opposite case of sort-range 
attractive surfaces, the adsorption is dominated by the 
short-range attraction. 

In Fig. 2a we present the phase diagram in the plane 
(csaitj c) differentiating the three regimes in the case of 
£ ps > 0, namely, when the surfactant-PE interaction is 
more favorable than the PE - bare surface interaction. 
For c < c*, in the area marked E, the short-range in- 
teraction between the surface and the PE chains is re- 
pulsive, and the adsorption is dominated by the electro- 
static attraction between the PEs and the surfactants. 
For high amounts of added salt, this electrostatic attrac- 
tion is screened, and the PE chains deplete (marked D). 
For c > c*, (marked C) the PE chains chemically adsorb 
on the surface for all c sa it. The electrostatic interactions 
in this region are dominated by the monomer- monomer 
electrostatic repulsion, and hence the addition of salt in- 
creases the adsorption, rather than causing depletion. In 
Fig. 2b we show the opposite case of e ps < 0. In this case 
the line separating the phases is no longer monotonic. 
For c < c*, the adsorption is mainly chemical, and is en- 
hanced by addition of salt. For c > c*, the adsorption 
becomes electrostatic and depends on c sa it. The amount 
of salt necessary for depletion decreases with c for low c 
since the increase in c increases the short-range repulsion 
between the surface and the PE chains, while the surface 
charge is not strong enough to electrostatically attract 
the PEs. For higher c, the increase in the surface charge 
with c increases mainly the electrostatic adsorption, and 
the necessary amount of salt increases with c. 



B. The Demixing Curve 



Repulsive Surface The spinodal lines of repulsive PE- 



surfaces (c* = 3.0, v ■ e ps = 0.56 2 a) for several values of 



Csait are presented in Fig. 3a. The spinodal of the pure 
surface monolayer, Eq. ©, is shown to be always higher 
than any of the charged-surface spinodals. All three 
spinodal lines are almost symmetric under the change 
c — > 1 — c, and the critical points are very close to the 
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FIG. 2: The numerically calculated phase diagram of the three PE- 
surfacant regimes are presented for two values of v ■ e pB and c*. In 
both figure parts the solid line differentiates the electrostatic adsorp- 
tion regime (marked E) from the depletion regime (D), and the dashed 
line differentiates both of them from the chemical adsorption regime 
(marked C). In (a) we set v ■ e ps = 0.5b 2 a and c* = 0.8, meaning 
that the short-range interaction of the PE chains with the charged 
surfactants is more favorable. For c < c* and low salt concentrations, 
the electrostatic attraction between the surfactants and the PE chains 
causes them to adsorb on the monolayer. For higher salt concentra- 
tions the PE and surface charges are screened by the salt ions, and 
the PE chains deplete. For c > c*, the short-range attraction between 
the surfactants and the PE chains is strong enough to allow PE ad- 
sorption for any salt concentration. In (b) we show the opposite case 
of v -£ ps = — 0.96 2 a and c* = 0.1, where the short-range interactions 
with the surfactants are less favorable. Here for c < c* the short- 
range interactions cause the PE to adsorb, while for c > c* the PE 
chains adsorb or deplete according to the strength of the electrostatic 
interactions and the added salt concentration. In both plots we use 
v = o.l. a = 5 A, v~ = 10 A 3 , <b1 = lO^A" 3 , Z = 1. f = 0.5. 
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FIG. 3: (a) The spinodal lines of a repulsive surface system are presented as function of the surface coverage c. The spinodal is shown for three 
amounts of added salt. The solid line corresponds to c sa i t = 0.1 M, the dashed line to c sa i t = 0.5 M and the dashed-dotted line to c sa i t = 1.0 M. 
All lines share c* = 3, a = 5 A, i; cx = 10 A 3 , <f>l = 10" 8 A" 3 , Z = 1, / = 1, 6 = 10 A, ue pB = 0.5b 2 a, T = 300 K, e = 80. All three spinodal 
lines are seen to be almost identical, differing only in the region close to the adsorption-depletion crossover, marked by an arrow. The dotted line 
marks the pure monolayer spinodal, Eq. and is shown for comparison, (b) The binodal (phase separation) line of a repulsive surface system 
is presented as function of the surface coverage c for c sa i t = 1.0 M. Other parameters are the same as in (a). The solid line is the binodal, while 
the dashed and dashed-dotted lines connect the separated phases for v = 0.2 and u = 0.31, respectively. The dotted line corresponds to the 
pure-monolayer binodal, Eq. J2J, and is shown for comparison. 



pure surface one at c c = 0.5. In Fig. 3b we present the 
binodal line of a repulsive surface for c sa i t = 1.0 M. The 
binodal line of the pure surface, Eq. is also presented, 
and is higher than the numerical binodal for all c. Like 
the spinodal line, the numerical binodal is almost sym- 
metric with respect to the change c — ► 1 — c. For high c 
values, all spinodal lines of Fig. 3a overlap almost iden- 
tically, showing that the dependence of the spinodal on 
the amount of added salt is very weak. For low c val- 
ues, careful observation of the three spinodal lines shows 
a difference between the lines, where the low salt lines 
are higher than the high salt ones. The three lines join 
together at c ~ 0.22 (see arrow). This difference be- 
tween the spinodal lines results from the PE depletion at 
low c. For low c values, the spinodal resembles a deple- 
tion spinodal, where addition of salt has a strong effect 
on the spinodal. For higher c values, the surface charge 
is high enough to allow PE adsorption. All three lines 
then move toward the electrostatic adsorption spinodal, 
that is weakly dependent on salt. The point where all 
lines join together is where the high salt spinodal crosses 
over from the depletion regime to the electrostatic ad- 
sorption regime, and its numerical value agrees with the 
adsorption-depletion crossover at 1.0 M of salt taken from 
Fig. 4. The spinodal lines in the different regimes are dis- 
cussed further in Sec. II VI 



In Fig. 4 the adsorbed amount T = 4% J °° da;(?7 2 — 1) 
is plotted as a function of the surfactant surface coverage 
c for twt = 1.0 M. We set v = 0.2, 0.31 in Figs. 4a, 4b, 
respectively. The adsorption is seen to be positive only 
from a threshold surface coverage Cd p , which is similar 
to the crossover point between spinodal lines shown in 
Fig. 3. Above this threshold, the adsorbed amount in- 
creases linearly with the surface coverage, in agreement 
with previous results [33, 34] . The threshold is found not 
to change with v. 

It is interesting to note the effect of the adsorption- 
depletion crossover on T in the coexistence region. In 
Fig. 3b we mark on the binodal line the coexisting phases 
for v = 0.2 (dashed line) and for v — 0.31 (dashed-dotted 
line). In Fig. 4 we show T for the same coexisting phases. 
Using the lever rule construction T inside the phase 
coexistence lies on a straight line between the two marked 
phases. For v = 0.2 (dashed line in Fig. 4a), the PE 
adsorbs in one of the phases while depleting from the 
other. As a result, the adsorption in the phase separated 
region is always higher than onto a reference single phase 
(solid line) , having the same mean surface coverage c, For 
v = 0.31 the PE adsorbs on the two coexisting phases. 
The adsorption in the phase separated region (dashed- 
dotted line in Fig. 4b) is almost identical to the adsorp- 
tion to the single corresponding mixed phase (solid line) , 
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and there is no gain in the adsorption due to demixing. 

The difference between these two behaviors can be 
qualitatively explained by the following argument. In 
the electrostatic adsorption case, the adsorbed amount 
of monomers scales linearly with the surface charge |30| . 
and hence also with c. The adsorption of the PE chains 
aims mainly to balance the surface charge, and since the 
demixing does not change the overall surface charge, the 
total adsorbed amount does not change as well. In the 
case of a depletion phase- adsorption phase coexistence, 
the charges on the low c (depletion) domain are screened 
by salt ions, while the higher surface charge domains are 
screened by the PE chains. In this case, even if the total 
surface charge is too weak to allow PE adsorption in a 
single mixed phase, the phase separation allows adsorp- 
tion at a part of the surface area, and thus allows for a 
positive adsorbed amount. The crossover between these 
two behaviors is at v ~ 0.3, corresponding to the low c 
phase being at c ~ Cd p . Namely, that the low c phase is 
at the adsorption-depletion transition. 

Short-Range Attractive Surface The opposite case of 
a short-range attractive surface is presented in Fig. 5. 
As can be seen from the figure, the spinodal line is much 
higher than the corresponding pure-surface spinodal line, 
and it is far less symmetrical. The addition of salt is 
shown to increase the spinodal temperature. This in- 
crease with salt can be explained by both the screen- 
ing of the electrostatic repulsion between the lipids, and 
by the increase in the adsorbed amount of monomers in 
this case (U , which increases the effective attraction be- 
tween the lipids. The corresponding adsorbed amount of 
monomers is shown in Fig. 6. The adsorbed amount here 
is convex in c, showing that the adsorption to a demixed 
two-phase system is always larger than to a single phase 
system. 

Attractive- Repulsive Crossover In Fig. 7 we show 
the spinodal lines of monomers for c* = 0.7 and both 
Csait = 1-0 M and c sa i t = 0.1 M. It is easily seen that 
the low salt spinodal is higher than the high salt one 
for c < c* (repulsive surface), while the opposite oc- 
curs for c > c*, consistent with Fig. 5. Although all 
three regimes occur in this spinodal, the curve has only 
one critical point, located near the crossover between the 
spinodal behavior of the electrostatic adsorption limit to 
the chemical adsorption limit, at c ~ c*. This shows that 
the transition of the spinodal line between the regimes 
is gradual rather than abrupt. The c sa i t = 1.0 M line 
changes convexity at c ~ 0.2, creating a "shoulder" in 
the spinodal. This shoulder marks the transition from 
a depletion regime to an electrostatic adsorption regime. 
For extremely strong short-range interactions (very high 
s ps ), this shoulder may transform into a triple point on 
the spinodal. However, we did not find such a triple point 
numerically for reasonable values of our parameter. 

In Fig. 8 the adsorbed amount for c* = 0.7 is plotted 
against the surface coverage c. For low c values, we see 
the PEs deplete, similar to Fig. 4, while for higher c the 
adsorption crosses over from a linear electrostatic adsorp- 



x 10" 




c\j 



0.2 0.4 0.6 0.8 1 



x 10 

■I i i i r 



(b) 




FIG. 4: The adsorbed amount V is plotted as a function of the 
surfactant coverage for two values of v. In (a) the solid line correspond 
toi/ = 0.2 while in (b) v = 0.31. Both figures share c salt = 1.0 M, 
and other parameters are the same as in Fig. 3. In both graphs, a 
low c depletion regime can be found, in which V is slightly negative. 
For larger surface coverage, the adsorbed amount scales linearly with 
the surface coverage, and no considerable convexity is obtained. The 
threshold surface coverage for adsorption is unaffected by v. The 
dashed line in (a) and dashed-dotted line in (b) are the adsorbed 
amounts inside the phase coexistence region, and are constructed by 
the lever rule. A significant increase in adsorption can only be achieved 
when one of the separated phases is a depletion phase. This is the 
case in (a). For higher values of v (b), the phase separated adsorption 
line is almost the same as the original (single phase) adsorption line, 
and no gain in adsorption is achieved. 
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FIG. 5: The spinodal line for a short-range attractive surface, and 
two salinities: c sa i t = 1 M (solid line) and c sa i t = 0.1 M (dashed line). 
The dotted line corresponds to the pure surface spinodal, Eq. 01 . 
The profiles were calculated for c* = —1 and Z = 0.1, and other 
parameters are the same as in Fig. 3. In contrast to the repulsive 
surface case, here the PE adsorption pushes the spinodal to higher 
temperatures than the spinodal of the pure monolayer. The effect of 
salt is also more pronounced than in the repulsive surface case, and 
it is shown that addition of salt increases the spinodal temperature 
considerably. 



FIG. 6: The adsorbed amounts of monomers are shown for the short- 
range attractive surface. The dashed line corresponds to c sa i t = 1 M, 
v = 0.3, the solid line to c sa lt = 0.5 M, v = 0.4 and the dashed-dotted 
line to c^t = 0.5 M, v = 0.5. All profiles share c* = 0, and other 
parameters are the same as in Figs. 3. In contrast to the repulsive 
surface case, here the addition of salt increases the adsorption. An 
increase in the temperature v decreases the adsorption as expected. 
The adsorption profiles T(c) are convex with respect to c, showing that 
a demixed lipid monolayer adsorbs more PEs than onto a homogeneous 
layer with the same mean coverage c. 



tion profile to a convex chemical adsorption profile. The 
addition of salt decreases the adsorption for low c val- 
ues, and increases it for higher surface coverage. Note, 
however, that the crossover between the two salt depen- 
dencies is not at c = c* but for a lower c. This result 
is in agreement with a previous one |3^ |. which shows 
that when the short-range repulsion between the surface 
and the PE chains becomes small enough, the adsorbed 
amount starts to increase with c sa it even before the sur- 
face becomes attractive, due to an increase in the surface 
monomer concentration rfs- 

So far we discussed the case where the short-range in- 
teractions between the PE chains and the surfactant are 
more favorable than with the bare surface (e ps > 0). In 
the case of a neutral polymer, changing the sign of e ps is 
exactly equivalent to replacing c* — * 1 — c* and c — > 1 — c. 
However, in our case the fact that the surfactant and 
polymers are oppositely charged breaks this symmetry. 

For very strong short-range interactions (large negative 
y-e vs ) and < c* < 0.5 , a triple point can emerge in the 
following way. For low c values, the spinodal resembles 
the attractive surface spinodal. The spinodal line in this 
region is always higher than the pure surface spinodal 
Eq. ©, as was shown in Fig. 5. For a P > c > c*, 
the spinodal line decreases substantially to the depletion 
spinodal, and for c > Cd p it moves to the electrostatic 
adsorption spinodal. Both of the latter spinodal lines 



are similar to those in Fig. 3a, both are lower than the 
pure surface spinodal, and have a critical point at c c ~ 
0.5. Therefore, there must be a triple point at c < 0.5, 
connecting three coexisting phases: (i) a low c chemical- 
adsorption phase; (ii) a low adsorption phase; and (iii) 
an electrostatic adsorption phase, as shown in Fig. 9. 

In Fig. 9 we show the numerically calculated triple 
point corresponding to the above conditions. At high 
temperatures v > 0.29, the phase coexistence for low 
surface coverage connects a chemical-adsorption phase 
with either another chemical adsorption phase, or a de- 
pletion phase. For higher surface coverage, the phase 
coexistence connects either a depletion phase or an elec- 
trostatic adsorption phase with another electrostatic ad- 
sorption phase. A triple point connecting all three phases 
occurs for c tr ~ 0.22, v iT — 0.29. 



IV. ANALYTICAL RESULTS 

In this section we present a theory to explain the three 
PE-surfactant regimes, and, in particular, aimed at ex- 
plaining the physics behind the spinodal behavior in each 
regime. We derive an approximate form for the free en- 
ergy in each regime, and show the resulting spinodal line. 
Our model is based on the scaling results of Refs. [30lf33| . 
and is shown to capture the correct behavior of the spin- 



8 



> 0.4- 




FIG. 7: The spinodal line is presented for c* = 0.7 for two salinities: 
Csalt = 0.1 M (dashed line) and c salt = 1.0 M (solid line). All other 
parameters are the same as in Fig. 3. The dotted line is the pure 
surface spinodal, and is presented for comparison. Both spinodals are 
very close, but one can see that the low salt spinodal line is slightly 
higher than the high salt spinodal line for c < c* , and slightly lower 
than it for c > c*. The saddle point on the high salt spinodal line 
corresponds to the adsorption-depletion crossover, and occurs at the 
same surface coverage as in the repulsive surface case. Both spinodals 
cross over from a low v range (lower than the pure surface spinodal) 
to a high v range (resembling an attractive surface) at c ~ 0.4. 
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x 10" 




FIG. 8: The adsorbed amount of monomers T is shown for c* = 0.7. 
The solid line corresponds to c sa i t = 0.5 M, v = 0.3, the dashed- 
dotted line to c aa j t = 1 M, v = 0.5 and the dashed line to c sa i t = 1 M, 
v = 0.3. All other parameters are the same as in Fig. 7. The increase 
in the amount of salt causes the adsorbed amount to decrease for low 
c, and to increase for high c. 



odal lines in each limit. The analytical equations can be 
compared to the numerical results with good accuracy, 
as shown below. 



A. Polyelectrolyte Depletion Regime 

When e ps (c — c*) < the surface repels the PE chains, 
and thus the only driving force for adsorption is the elec- 
trostatic attraction between the surfactants and the PE 
chains. This electrostatic attraction strongly depends on 
the amount of small ions in the solution [23] ■ When this 
amount increases beyond a threshold value, the salt ions 
neutralize the surface and PE charges, and the solution 
behaves like a neutral polymer solution. In this case, the 
repulsive short-range interaction between the PE chains 
and the surface causes the PE chains to desorb entirely 
from the surface. This threshold value of salt was previ- 
ously shown [33| to depend on the surface charge density 
of the interface and on the charge of the PEs via the scal- 



-V2, 



This relation can be 



ing law (c s * alt ) 3 / 2 ~ \a\ fl 
inverted to give a threshold surfactant concentration Cd p : 



Cdp 



A 



3/2^1/2^2^2 



dp" 



'-salt'S 



Zf 



(12) 



below which the PE chains do not adsorb on the charged 
surface. The prefactor Ad P can be calculated by compar- 
ison of Eq. 1|12[) with the numerical results as in Figs. 3-4, 
yielding Ad p — 2.35. 

When the PE chains desorb from the surface, the elec- 
trostatic potential depends only on the adsorption profile 
of the small ions. According to the well known Poisson- 
Boltzmann theory |40j , this decay is exponential as func- 
tion of the distance from the surface x. This exponential 
decay causes the interaction term /<^f?j 2 C between the 
PE chains and the electrostatic potential in Eq. © to be 
small in comparison to the other terms. The free energy 
hence decouples and describes a neutral polymer and an 
ionic solution. Using the Debye-Hiickel approximation 
for the small electrostatic potentials, the total free en- 
ergy becomes: 



F 



Kurf 



Kb 4 



(13) 



a 2 v ex c 



18e ps (c* - c) 



b 4 a 2 v cx (f) 2 
3£p S (c* - c) 



The second term is the free energy of a solution of small 
ions under the Debye-Hiickel approximation, the third 
term is the free energy of a neutral polymer solution as 
in Ref. 

the neutral surfactant monolayer 



and the first term -F s ' urf is the free energy of 



F s ' urf = b~ 2 [cine + (1 - c) In (1 - c) 



(14) 



+u~ 1 c(l - c) 
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where we note that the surface charge term was already 
taken into account in the small ion solution term of 
Eq. 11311 . Differentiating this free energy twice with re- 
spect to c, we obtain the following equation for the spin- 
odal: 



1 -c 



2vl 



18e ps (c* - c) 



AttI b Z 2 

+ Kb 4 

£p S (c* - c) 2 



(15) 



The first term of the spinodal Eq. l(T5|l is the pure surfac- 
tant monolayer spinodal as in Eq. J2J), while the second 
term accounts for the electrostatic interactions between 
the salt ions and the surface and the third term for the 
neutral polymer solution 3] . The neutral polymer term, 
as was shown in Ref. Q, causes the spinodal line to in- 
crease slightly. The magnitude of this term is inversely 
proportional to the short-range interaction strength £ ps , 
and can be discarded in the limit of strong short-range re- 
pulsion. The salt term, in turn, pushes the spinodal line 
downwards, as was previously shown by Ref. 4, 5]. The 
decrease in the spinodal line is substantial, and depends 
on the amount of added salt. In this (Debye Hiickel) 
approximation, for very large amounts of added salt all 
electrostatic interactions between the monomers and the 
surfactants are totally screened, resulting in a strong de- 
crease in F c \. The spinodal line then fully recovers the 
neutral polymer limit from Ref. Q. For lower amounts 
of salt, which include most physical cases, the spinodal 
is pushed downwards from the pure monolayer spinodal, 
as seen in Fig. 3a. 




0.2 0.4 0.6 0.8 1 



FIG. 9: Binodal and spinodal lines exhibiting a triple point. The 
dashed-dotted line is a numerically calculated spinodal line, for the 



following parameters: a = 5 A, b = 
Z = 1, c* = 0.1 = 1 M, T 



10 A, V ■ £ D 



-1.06 2 a, / = 0.3, 
10~ 8 A~ 3 , d cx = 



10 A . The solid line represents the corresponding binodal line. A 
triple point, connecting three coexisting phases, exists for utr — 0.29. 
The three coexisting phases are connected by the dashed line. 



width (denoted by D). The free energy of the solution 
can be calculated using a Flory-like approximation [30| . 
While the full details can be found in the previous pub- 
lications, we present here the main results: 



B. Electrostatic Adsorption Regime 



For c > Cd p and s ps (c — c*) < 0, the PEs are attracted 
to the surface only via their electrostatic interactions, 
but cannot approach the surface because of the short- 
range repulsive interactions. Due to the exclusion from 
the surface, the PE surface concentration becomes very 
small, and can be ignored at first a ppr oximation. 

In previous publications [30L l33l |34 | . we discussed at 
length the adsorption of polyelectrolytes to charged sur- 
faces, which have short-range repulsive interactions with 
the PE chains. We summarize here the main results. It 
was shown that the dominant interaction is the electro- 
static interaction between the PE and the surface. The 
excluded volume interaction was found to be negligible in 
the numerical solutions of the mean-field Eqs. (|8I9|I . The 
PE adsorption was shown to decrease with the addition 
of salt [33, y4| for low amounts of added salt, while for 
higher amounts the polymer chains were shown to des- 
orb entirely from the surface |3^| . It was shown that the 
adsorbed polyelectrolyte layer can be characterized by a 
single (rescaled) monomer concentration scale (denoted 
by r/Jj) and a single length scale for the adsorbed layer 



g| = Al "_^L A 2f | ff | I bD ^Wm + (16) 
AsflB^tvifD 3 + A 4 \a\ 2 l B D 

where the surface charge is denoted by a — —cZ/b 2 . The 
first term is the elasticity term of the polyelectrolyte, the 
second is the electrostatic attraction between the surface 
and the monomers, and the third is the electrostatic re- 
pulsion between monomers. This repulsion is assumed 
to be strong enough to dominate over the excluded vol- 
ume interactions. The last term is the interaction of the 
electrostatic potential and the surface charge, where the 
rescaled potential is |£ s | = A4 \a\ IbD. The short-range 
interaction between the surfactants and the PE chains de- 
pends on the surface monomer concentration, and can be 
ignored because of the strong short-range repulsion. Min- 
imization of Eq. (|16fl with respect to D and r^ M yields: 
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D 



J.2 2 



1/3 



Infer 
a 2 Vh 



af 



2/3 



-Fads — Iq 



2,2 \ 1/3 



a 2 1 



f 



,5/3 



(17) 
(18) 
(19) 



where the prefactor Iq can be shown to be positive. In 
this case, the full free energy can be written as: 



F = In 



a 2 l 2 B Z 5 

10 



fb 



1/3 



„5/3 



(20) 



cLr±c+ (1 - c) ln(l - c) + v~ x c(\ - c) 



Note that this free energy has a fractional power (non- 
analytic) in c. The spinodal equation is obtained by dif- 
ferentiating Eq. (|20[1 twice with respect to c, yielding: 



same as that of the monomer concentration. The single 
length scale allows us to use similar scaling ideas as in 
the repulsive surface case. The highest monomer concen- 
tration occurs at the surface itself 77 ma x = Vs, and acts 
as the scale for the monomer concentration. Similarly to 
Ref. [3(j, we can produce scaling rules for the monomer 
adsorption using a Flory-like free energy. The free energy 
in this case is: 



-Fads 



2^2 



A 1 



6 D 



-ps ( c 

-A 3 f 2 h 



(23) 



where the first term is again the contribution of the chain 
elasticity, the second term is the short-range attraction 
of the polymer chains to the surface, and the last term is 
the electrostatic repulsion between the monomer layers. 
The electrostatic interactions between the surface and 
the monomer concentrations are assumed to be small in 
comparison to the short-range ones, and are omitted from 
this free energy. Minimization of this free energy yields 
the following scales for D and rj 2 : 



2^ = - 
v c 



1 - c 



10 J 



a 2;| z 5 X l/3 



fb 4 



.-1/3 



(21) 



Eq. Q21JI shows that the spinodal temperature v sp is al- 
ways lower than that of the pure surfactant monolayer, 
Eq. (PJ). This result is consistent with results derived for 
the spinodals of protein solutions The critical point 

is the maximum of this spinodal line f sp (c) in Eq. (|21|l . 
satisfying: 



2c, 



cl (1 



-27 J ° 



a 2 l 2 B Z 5 
fb 4 



1/3 



-4/3 



(22) 



The critical point is at c > 0.5, regardless of the sign of 
e ps . For our chosen parameters of a = 5 A, Iq = 2.8, 
Ib — 7 A, b — 10 A and Z = 1 the critical concentration 
from Eq. l(2*2^l turns out to be at c = 0.5365, showing 
that the spinodal is almost symmetric, in agreement with 
Fig. 3. 



C. Strong Short-Range Attraction Regime 

For large enough e ps (c — c*) > (strong short range 
attraction) the adsorbing polymer layer is highly charged, 
more than the initial surface [34l l35| . In this case, the 
main electrostatic contribution to the free energy comes 
from the monomer-monomer repulsion, rather than the 
electrostatic attraction to the surface. 

When the amount of small ions in the solution is small, 
the length scale for electrostatic potential decay is the 



2^2 



D 

^ads 

kB~f 



a 2 b 



£ps (C-C*) 



£ ps (C 



b s aH B f 2 ^ 2 



b w aH B p 



(24) 
(25) 
(26) 



where I\ is a constant of order unity. The adsorption free 
energy is negative, as expected, and is very strongly de- 
pendent on the strength of the short range interactions. 
The total free energy is now the sum of the pure mono- 
layer free energy and the adsorption free energy, and the 
spinodal line can be readily obtained as: 



20L 



9 5 (c-c*)- 
^ p l B a%8f 



(27) 



where we denote = v ■ e ps . In this case, the spinodal 
temperature v> sp is always larger than the corresponding 
temperature of the pure-surface spinodal. The critical 
point is obtained by differentiating Eq. (|27J) once more 
with respect to c. 



(1-Cc)' 



-60J r 



e 5 



c *y 



vll B a%*f 2 



= 



(28) 



which for 9 > gives c c > 0.5, and vice versa. The 
resulting spinodal line is less symmetric than in the 
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FIG. 10: Comparison between the numerical and analytical spinodal lines. In (a) we present the numerical spinodal lines for two sets of 
parameters. The solid line is the numerical spinodal line in the repulsive case for Z = 1, c* = 3.0, while the dashed line is the numerical spinodal 
line in the attractive case, for Z = 0.1, c* = —1.0. Both spinodal lines share a = 5A, 6 = 10A, v ■ e ps = 0.5b 2 a, / = 1, c sa i t = 0.01 M, 
T = 300 K, 0jj = 10 — 8 A - 3 , t; cx = 10 A 3 . We do not show the analytical spinodals here since they are too close to the numerical ones to be 
drawn. In (b) we present the differences between the analytical calculations and the numerical spinodals. The solid and dashed lines are the 
difference between the spinodal lines in Eqs. 1211 and 1271 . respectively, with Iq = 2.8 and I\ = 0.023, and the numerical spinodals in (a). The 
differences are small, and the numerical spinodal lines are well characterized by the analytical approximation. 



electrostatic adsorption case, and the critical point is fur- 
ther from the pure surface value of 0.5 than the repulsive- 
surface spinodal Eq i|21|) ■ This result is also in agreement 
with the numerical results of Fig. 5. 



D. Comparison of the analytical and numerical 
results 

Both spinodals, Eqs. (|21l27p . contain one free param- 
eter each, and can be compared to the numerical spin- 
odals. The comparison is presented in Fig. 10. As can 
be seen, both of the spinodal equations show good agree- 
ment with the numerical lines, using the values Iq = 2.8 
for the repulsive surface case and I\ = 0.023 for the at- 
tractive surface case. 

The analytical model given above for the two adsorp- 
tion regimes is accurate only for the case of low amounts 
of added salt. This is because the main assumption is 
the existence of a single length scale for the decay of the 
electrostatic potential and the monomer concentration. 
More work is necessary to extend the present model for 
the high salt regime. 



V. CONCLUSIONS 

We present numerical and analytical results for the 
phase separation of a surfactant monolayer in presence 
of polyelectrolytes (PEs). We show that the short-range 
interactions between the surfactant and polyelectrolytes 
have a strong impact on the surfactant critical tempera- 
ture and critical concentration. We also present analyti- 
cal results in two limits: (i) strong short-range repulsive, 
and (ii) attractive surfactant-polyelectrolyte interactions. 

For short-range repulsive interactions, our results show 
that surfactant phase separation occurs for lower temper- 
atures than for a non-interacting surfactant layer. For 
such systems, the polyelectrolyte chains may desorb com- 
pletely from the low surfactant concentration domains, 
while for higher concentration phases the PE adsorp- 
tion increase linearly with the surfactant surface cover- 
age. The critical surfactant concentration slightly devi- 
ates from the bare one c c = 0.5, showing nearly symmet- 
rical spinodal and binodal lines. In addition, we hnd a 
temperature range in which phase separation does not 
increase the PE adsorption. For lower temperatures, 
where the phase separation occurs between a strongly 
adsorbing phase and a depleting phase, the gain in ad- 
sorption becomes substantial. For short-range attrac- 
tive interactions, in contrast, the phase separation al- 
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ways increases adsorption, and the demixing tempera- 
tures are always larger than the pure-surface ones. The 
critical point moves more significantly away from the 
pure-surface symmetric value. The convexity of the ad- 
sorption is more pronounced, showing that phase separa- 
tion indeed increases the PE adsorption on the surfactant 
monolayer. In intermediate cases, where the short-range 
interaction between the polyclectrolyte chains and the 
surface changes sign as function of the surfactant surface 
coverages, a triple point can be found on the binodal line 
connecting three coexisting phases on the same surface. 

The addition of salt is shown to increase the threshold 
surfactant-concentration for the desorption phase. It has 
a weaker effect on the spinodal line for the case of elec- 
trostatic adsorption. For short-range adsorbing surfaces, 
the addition of salt increases the spinodal line consider- 
ably. 

Our results are derived within mean-field, using the 
ground state dominance approximation. This treatment 



neglects correlation effects between the different con- 
stituents, and they may be strong in some cases. We 
also assume that the lipids are entirely insoluble, and 
neglect the effects of the line tension between the surfac- 
tant phases. Despite these drawbacks, we believe that 
our main results still hold. Our results may serve as a 
starting point for more complex models. 
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